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Abstract 



o 

In an extended, new form of general relativity, which is a teleparallel theory 

O of gravity, we examine the energy-momentum and angular momentum carried 

O" 1 

^ ■ by gravitational wave radiated from Newtonian point masses in a weak-field 

approximation. The resulting wave form is identical to the corresponding 
wave form in general relativity, which is consistent with previous results in 
teleparallel theory. The expression for the dynamical energy-momentum den- 
sity is identical to that for the canonical energy-momentum density in general 
relativity up to leading order terms on the boundary of a large sphere includ- 
ing the gravitational source, and the loss of dynamical energy-momentum, 
which is the generator of internal translations, is the same as that of the 
canonical energy-momentum in general relativity. Under certain asymptotic 
conditions for a non-dynamical Higgs-type field tp k , the loss of "spin" angular 
momentum, which is the generator of internal SL(2, C) transformations, is 
the same as that of angular momentum in general relativity, and the losses 
of canonical energy-momentum and orbital angular momentum, which con- 
stitute the generator of Poincare coordinate transformations, are vanishing. 
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The results indicate that our definitions of the dynamical energy-momentum 
and angular momentum densities in this extended new general relativity work 
well for gravitational wave radiations, and the extended new general relativity 
accounts for the Hulse- Taylor measurement of the pulsar PSR1913+16. 



1 Introduction 

General relativity (GR) is a standard theory of gravity which has passed all the 
observational tests so far carried out, and it constitutes, together with quantum 
field theory, a basic framework of modern theoretical physics. In GR, however, it 
is usually asserted [|l| that well-behaved energy-momentum and angular momentum 
densities cannot be defined in general for a gravitational field. For a restricted class of 
systems including asymptotically flat space-time, there exist tensor densities whose 
integrals over the cross section of the null infinity give the energy-momentum and 
angular momentum of the system in question ||. 

There are many theories Pfl that are potential alternatives to GR, including the 



Poincare gauge theory (PGT) Q and extended new general relativity (ENGR) [§. 
PGT is formulated on the basis of the principal fiber bundle over the space-time 
possessing the covering group Pq of the proper orthochronous Poincare group as the 
structure group, following the standard geometric formulation of Yang-Mills theories 
as closely as possible. ENGR is formulated as the teleparallel limit of PGT. The 
dynamical energy-momentum and "spin" angular momentum densities of gravita- 
tional and matter fields are all space-time vector densities, and their integrations 
over an arbitrary space-like surface a are well-defined for any coordinate system 
employed ||. 

For asymptotically flat space-time whose vierbeins satisfy certain asymptotic 
conditions,^ the integration of the dynamical energy-momentum density over o is 
the generator of internal translations and gives the total energy-momentum of the 
system. Also, the integration of the "spin" angular momentum density over a is the 
generator of internal SL(2, (^-transformations and gives the total (=spin+ orbital) 



J With regard to the situation in ENGR, see Eqs. ( |2.45D and fl2.46| ) in §U 
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angular momentum in both theories. This holds in PGT when the Higgs-type field 
tjj k satisfies the asymptotic condition ip k = e^^x^ + ip^ h + 0(l/r@) with con- 
stants e^ k ^, ip^ k J7|, || ^] and in ENGR when this asymptotic condition and certain 
other additional conditions are satisfied 0. These theories describe within the un- 
certainties all the observed gravitational phenomena when the parameters in the 
gravitational Lagrangian densities satisfy certain conditions. 

Direct observation of gravitational waves is one of the most challenging problems 
in present day gravitational physics. Several projects designed for this purpose are 
now being carried out, and gravitational radiations from various possible sources 
have been investigated theoretically, mainly on the basis of GR. However, also in 
classes of teleparallel theories of gravity, the form of gravitational waves is known 



to be identical to that of GR in post-Newtonian approximations [10, 11 



For the case in which a gravitational wave is radiated, however, the asymptotic 
behavior of vierbeins is different from that considered in Ref. || in general, and the 
question of whether our definitions of the energy-momentum and angular momentum 
densities work well in this case should also be answered. 

The purpose of the present paper is to examine, in a weak-field approximation, 
the energy-momentum and angular momentum carried by gravitational waves radi- 
ated from Newtonian point masses. In §0, the basic framework of ENGR is briefly 
summarized as preparation for later discussion. In the forms of the gravitational 
field equations and the dynamical energy-momentum density T k ^ of the gravita- 
tional field are given in the weak-field approximation. For plane wave solutions of 
the linearized homogeneous equations of a gravitational field, we give the average 
of G T k fl over a space-time region much larger than the inverse of the absolute value 
of the three-dimensional wave number vector. In §§, the quadrupole radiation for- 
mula for a gravitational wave emitted from a system of Newtonian point masses is 
obtained. In §|5], we examine the emission rates of the dynamical energy-momentum 
and the angular momentum for two types of the asymptotic form of a Higgs-type 
field. Further, the emission rates of the canonical energy-momentum and the "ex- 
tended orbital angular momentum" are examined. Finally, in §|| we give a summary 
and discussion. 
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2 Basic framework of extended new general rela- 
tivity 



2.1 Poincare gauge theory 

We first give the outline of PGT, because ENGR is formulated as a reduction of 
this theory. 

PGT is formulated on the basis of the principal fiber bundle V over the space- 
time M possessing the covering group Pq of the proper orthochronous Poincare group 
as the structure group. The space-time M is assumed to be a noncompact four- 
dimensional different iable manifold with a countable base. The bundle V admits 
a connection T, the translational and rotational parts of whose coefficients will be 
written A k [1 and A h lfl , respectively. The fundamental field variables are A k [i and A h lfl , 
the Higgs-type field is ip = {i> k }, and the matter field is = — 1,2, ... , iV}|] 

These fields transform according tof] 

^^{Kia-^f^-n, (2.1a) 

A'\ = (A(a- 1 )) fe m (A™ + f™ + A m n /), (2.1b) 

A'% = (A(a- 1 )) fe m A m „ M (A(a))^ + (A(a- 1 )) fc m (A(a))^ , (2.1c) 

<P' A = [p((t,a)- 1 )] A B( p B , (2.1d) 



under the Poincare gauge transformation 

cr'(x) = a{x) ■ [t(x),a(x)] , 

t(x)eT\ a(x) G SL(2, C). (2.2) 

Here, A is the covering map from SL(2, C) to the proper orthochronous Lorentz 

group, and p denotes the representation of the Poincare group to which the field <fr 

5 Unless otherwise stated, we use the following conventions for indices. Letters from the middle 
part of the Greek alphabet, A, fj,, . . ., and from the middle part of the Latin alphabet, k, I, m, 
. . ., take the values 0, 1, 2 and 3. The capital letters A and B are used as indices for components 
of the field <f>, and N denotes the dimension of the representation p. 



6 For the function / on M, we define / jM = df/dx 11 . 
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belongs. Also, o and a' stand for local cross sections of V. The dual components e k fl 
of the vierbein fields e^d/dx^ are related to the field ip k and the gauge potentials 
A k ^ and A k ^ through the relation 

e\ = ^ + A%rlt + A\ , (2.3) 

and these transform according to 

e% = (ACa- 1 ))*^ , (2.4) 

under the transformation (|2.2|). Also, they are related to the metric g^dx^ ® dx u 
of M through the relation 

g^u = e\r) M e l u , (2.5) 

with (r] kl ) = diag(-l, 1, 1, 1). 

The field strengths R k lfiu , R k flL , and T k of A k lfl , A k and e k fl are given byQ 

I&inv = ^(A k i[ u ^] + A k m ^A m lu ^), (2.6a) 
R k lxu d ^2(A k M +A k l[fX A l v] ), (2.6b) 

^ = 2(^+^1)- ( 2 - 6c ) 
respectively, and we have the relation 

T\ v = R k ^ + R\^ 1 . (2.7) 

The field strengths T ^ and R \ are both invariant under internal translations. 

There is a 2 to 1 bundle homomorphism F from V to the affine frame bundle 
A(M) over M, and there exist an extended spinor structure and a spinor structure 



associated with it [12]. The space-time M is orientable, which follows from its 



assumed noncompactness and the fact that M has a spinor structure. 



'We define 

dof 

2 



A. ..[fj,.. .„].., — -(A../i...i/... — A. .v... I*...), 
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The affine frame bundle A(M) admits a connection Fa- The T 4 part and the 
GL(4, R) part Y x fiu of its connection coefficients are related to A k lfl and e k fl through 
the relations 

= 5*, , (2.8a) 
A% = e\e\T\^ + e\e\^ , (2.8b) 

by the requirement that F maps the connection T into Fa, and the space-time M 
is of the Riemann-Cartan type. 
The torsion is given by 

T\ v ^ 2r A M , (2.9) 
and the T 4 and GL(4, R) parts of the curvature are given by 

R x =2{Y x +Y x Y p v X (2.10) 



pv V [u,n] ' 1 p[p L u]J 

= 2(r x 4- r A r r 

ppv \ p[v,p] T \P P v \' 



R X puv — 2(r A r i + r A T r r r ,), (2.11, 



respectively. Then, we have the relations 



= e\T\ v = e\R x , (2.12) 



pv c A fiv c A L /i^ ) 
_ p fc p P E>A 

i/Lty c r l Jl ppv i 



R\»u = e\e P iR* , (2-13) 



which follow from Eq. 

The covariant derivative of the matter field takes the form 

D k <j) A = e\D^ A , (2.14) 
D,<p A = d,<p A + l -A kl ^{M kl 4>) A + iA\{P k <t>) A , (2.15) 

where M k i and P k are representation matrices of the standard basis of the Lie algebra 
of the group Pq\ Mki = —ip*(M k i), Pk = —ip*(Pk)- The matrix P k represents the 
"intrinsic energy-momentum" of the field (p A ]12| , and it is vanishing for all observed 
fields. 
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2.2 Extended new general relativity 

In PGT, we consider the case in which the field strength R kl ^ vanishes identically, 

R k \ v = 0. (2.16) 



Thus, the curvature R x pul/ vanishes, and we have a teleparallel theory. 



p\i,v 

By choosing the SL(2, C)-gauge such that 

A k \ ee 0, (2.17) 
the following reduced expressions are obtained: 

e\ = i)\ + A\, (2.18) 

rV = > ( 2 - 19 ) 

D k <j) A = e\D^ A , (2.20) 
D^ A = d^ A + iA\{P h <t>) A . (2.21) 

The Lagrangian takes the formQ 

L = L T {T klm ) + L M (e\, ^, D k( p A , A ), (2.22) 

where L M is the Lagrangian of the matter field <p A and L T is the gravitational 
Lagrangian. We impose the following requirements: (R.i) L is invariant under the 
transformation ( |2.2| ) with arbitrary functions t k and an arbitrary constant element 
a of SL(2, C); (R.ii) The functional L is a scalar field on M. 
The gravitational Lagrangian [I3| 

L T dcf Clt Mm tHm + + (2>23) 

satisfies these requirements, where Ci, c 2 and c 3 are real constants. The quantities 
tkim, v k and a k are the irreducible components of T k i m defined by 

f 1 1 1 

tkim = -(T k lm + Tikm) + -(rjmkVi + r] ml V k ) - -TjuVm , (2.24) 

2d 6 
v k = T\ k , (2.25) 

Ofc = j;CkimnT lmn , (2.26) 
6 



8 The field components e fc M and e M fc are used to convert Latin and Greek indices. Also, the raising 
and lowering of the indices fc, I, m, . . . are accomplished through use of (rj kl ) = (ffez) -1 and 
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where the symbol eki mn represents the Levi-Civita tensor, with e(o)(i)(2)(3) = — 1-Q 
If the parameters ci, C2 and C3 satisfy 

ci = -02 = ^c 3 = , (2.27) 
where k is the Einstein gravitational constant, k == 87rG/c 4 J^ the gravitational part 



of the action integral is equal to the Einstein-Hilbert action integral, namely [JIB 



J d 4 xV^L T = J d 4 x±^R({}), (2.28) 



where R{{}) denotes the Riemann-Christoffel scalar curvature. Here, we should 
mention that even in the case that the condition ( |2.27| ) is satisfied, our theory 
does not reduce to GR, because the couplings of matter fields (the spinor field, for 
example) with the gravitational field are different from those in GR. 

In defining the energy-momentum and angular momentum, there are two pos- 
sibilities for choosing the set of independent field variables H 0, [|: O ne is to 
choose the set {ip k , A k ^, 4> A }, and the other is to choose the set {ip k , G k , 4> A }- In the 
rest of this paper, we employ {ip k , A k ^, <f) A } as the set of independent field variables, 
because this choice is superior to the other, as shown in Refs. || [7], j|, || . 

From the requirement (R.i), we obtain the identities^ 

w + ^& +i w^ )As0 - (2 ' 29) 

F k ifiu) = 0, (2.30) 
to V ~ W ~ jjjr = 0, (2-31) 

d ^ s » - ~ 2 w; Al] » _ l w [Mkl(p)A " °' (2 - 32) 



9 Latin indices are put in parentheses to distinguish them from Greek indices. 
10 G and c stand for the Newtonian gravitational constant and the light velocity in vacuum, 
respectively. 

n For instance, 5L/Sip k denotes the Euler derivative with respect to ip k . 
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where we have defined 



T < ~w; +i m (m) - (2 - 35) 



L = V^gL , g = det(^) , (2.33) 
= ^ • (2-34) 

"""" 9L - +i w} mA 

tot S k f = - 2 VAl, - ^ (Mfe ^ )A • (2 ' 36) 

If the field equations 6L/6A k fl = and 5L/5(p A = are both satisfied, we have 
the following: 

• The field equation 5L/5ip k = is automatically satisfied, and thus ip k is n °t 
an independent dynamical field variable. 

• There are two conservation laws, 

d^Tjf = 0, (2.37) 
V* V = 0, (2.38) 
which follow from Eqs. (g^)-(ggg). 

The former is the differential conservation law of the dynamical energy- 
momentum, and the latter is that of the "spin" angular momentum. 

We split the densities tot T fc M and tot ^ M into gravitational and matter parts as 

tot S k f = G S k f + M S k f , (2.40) 

where we have defined 

°V = - 2Jl/%, , (2.43) 
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with L T = f y/—gL T and L M = f ^J—gL M . Here, G T k M and M T, M are the dynamical 



- fe -*■ k 

energy-momentum densities of the gravitational field and the matter field, respec- 
tively, while G S kl fJ ' and M are the "spin" angular momentum densities of the 
gravitational field and the matter field, respectively. The densities T^, M T fc M , 
G S kl tM and M are all space-time vector densities 0. 

In Ref . || , the integrals of the dynamical energy-momentum and "spin" angular 
momentum densities over a space-like surface a are examined for vierbeins with the 
asymptotic behavior described below. 

(1) The components e k /1 of the vierbein fields possess the asymptotic property^ 

e\ = e^\ + f k ^ f\ {m) = 0{l/r l+m ) (m = 0,l,2), (2.45) 



where jf^( m ) denotes the mth order partial derivative with respect to x x , and the 

e^ k ^ are constant vierbeins satisfying e^ k ^r] k [e^ 1 v = rj^. 

(2) The antisymmetric part of components f^ v = f e iflkif 1 !, satisfy 



/[H,(m) = 0(l/r 1+a+m ) (m = 0, 1), (2.46) 



where a is positive but otherwise arbitrary. 
It has been shown that 



M k = I tot T k »da, = e {Q \M, , (2.47) 

J a 

Sh = [ tt V^ = e ( V (0) ^ + W ( \ M i] > ( 2 - 48 ) 
where da^ denotes the surface element on a. In the above, we have defined 

M M d ^ f Vflv ! 9 vX da x (2.49) 

J a 

= ( d p K^ xp da x = j {x»6 uX - x u 6» x )d<j x , (2.50) 



12 The expression 0(l/r n ) with real n denotes a term for which r n O(l/r n ) remains finite for 
r — > oo; a term denoted as 0(l/r n ) could, of course, also be zero. 
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with 



it/A def 1 



d p d a {(-g)g^g^} 



(2.51) 



K^p d ^ f I (^{(- g )^[V la } " ^<9 <T {(-^ /4 V ,T } + (-ff)^ [ V ] " - rf [X rf ]v ) ■ 

(2.52) 

This expression of 9 uX is the same as that of the symmetric energy-momentum 
density proposed by Landau and Lifshitz [13]. Equation ( |2.48| ) has been obtained 
by choosing the asymptotic form of the Higgs-type field ip k as 

(2.53a) 
(2.53b) 
(2.53c) 



^ = e (°) fe ^ + 4>(° )k + o(-^), 



T (3 

09 >0) 



k = o 



with and (3 constant, whereas Eq. ( 2.47 ) has been obtained without imposing 
the conditions ( [2.53j ). 

From the requirement (R.ii), we obtain the identity 



with 



™ v def , VT jp Xu A k 9L A dL 



i^A def j-, uXj^k 



A;/ 



The identity fl2.54| ) leads to 



3 v t; = o, 

d x M u vX = 



when 5L/5A k u = 0, where we have defined 



M 



v\ def 



2(*; A -x-z; A ). 
n 



(2.54) 

(2.55) 
(2.56) 



(2.57) 
(2.58) 



(2.59) 



Equations ( 2.57|) and ( 2.58 ) are the differential conservation laws of the canonical 
energy-momentum and the "extended orbital angular momentum" defined by 

M C M ^ I f;da u , L; d = f I M; x da x , (2.60) 
J & J (j 

respectively . The canonical energy- momentum and the "extended orbital angular 
momentum" are the generators of general affine coordinate transformations. The 
antisymmetric part == L^1]\ U ] is the orbital angular momentump 3 ] and is the 
generator of coordinate Lorentz transformation ||. 

We split the canonical energy-momentum density into gravitational and matter 
parts as 

f/ = G f/ + M T/, (2.61) 

where we have defined 



G T u def 5;L T _ ^ _ ^ k ^ (2 _g 2) 



f)T M f)T M 

=5 » L -W?»~W?"' { ] 

The density G T p u does not transform as a tensor density under general coordinate 
transformations, while M T p does transform as a tensor density. Q 

As is described in Ref. ||, the generators M c and L " vanish for vierbeins with 
the asymptotic forms satisfying Eqs. ( |2.45| ) and (|2.46|) when the condition ( |2.53|) is 
satisfied. 

The field equation 5L/5A k ll = has the expression 

-2V A i^ A + 2v x F» uX + 2H» V - g pu L T = , (2.64) 

where we have defined 

V A F^ A = d x F^ x + r» aX F™ x + T u aX F^ x + T\ x F^ a , (2.65) 

1 
3 



F „u\ drf Ci ^ux _ t ^ + C2 f g ^ v X _ g ^ v ^ _ I C3 e^Ap ap ; (2.66) 



H*" = T p(Jp F prj v - -T^F^ . (2.67) 



13 This Ir^j should not be confused with the orbital part of the "spin" angular momentum Ski- 
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Also, T^ u is the energy-momentum density of the gravitational source defined by 

5L M 
~5AK ' 



XT M 

' = WiThr ■ ( 2 - 68 ) 



3 Weak-field approximation 



We now consider weak field situations in which the vierbein fields e k ,, take the form 



e 



k _ J0)k , fk I fk 



; + i/*j<i, (3.i) 



where the e^ ""^ are constant vierbeins satisfying e^ k ^r]kie^ 1 u = f]^. The compo- 
nents of the metric and torsion tensors are given by, up to terms linear in /* 

= rip, + 2f^) , (3.2) 
T\fj,u = dfj,f\ u — d u fx[i , (3.3) 

where we have defined f^ v = f e^^jw/'i/Q 

In the weak-field approximation, the symmetric and antisymmetric parts of the 
field equation ( |2.64| ) take the form 

3d {Df ipu) - d\dj {v x) + dJw) + V^d P dJ^} 

+ (ci + c 2 ) {-v^nf - 2 % „d p dJ^ + d,dj 

+ d\dJ {vX) + dj^x)) - d\d,f wx] + dj^x])} = T {fiv) , (3.4) 



9 

+ (ci + c 2 ) {d\dJ {uX ) ~ dJw) - d\d,f W x] - d„f w )} = T [H , (3.5) 



with □ = f d^d^. Here, we have introduced 



/(H = f iW) - , f = yTfu*) , (3-6) 

/ = • (3.7) 
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For /* we use the convention that both the Greek and Latin indices of /* are raised or 
lowered with the Minkowski metric, and that they are converted into one another with e'' ^ or 
e (0)/ ;, where (e (0)/ £) d = (e (0) * Thus, for example, represents r) vX e {0 \f\& g vX e\f\). 



13 



We consider the energy-momentum density of the source to lowest order in 
f k ^. Therefore it is independent of f k and satisfies the ordinary conservation law in 
special relativity, 



d v T, 



0. 



Let us consider the transformations 



f(ljLV) ~ f(^) d{i. E V dyE^ , 



(3.8) 

(3.9) 
(3.10) 



where and Xn are arbitrary small functions. Since Eqs. ( |3.4| ) and ( ^.5| ) are in- 
variant under the transformations (|3.9|) and (|3.10|) with = x M , we can impose the 
harmonic coordinate condition 



v /CH = 0. 

The Lagrangian L T with the parameters c\ and C2 satisfying 



Cl 



-c 2 = -T- 
3k 



(3.11) 



(3.12) 



compares quite favorably with experiment ||13|| . We therefore assume ( |3.12j ) to hold 
henceforth. Under the conditions ( |3.11| ) and ( |3.12|) , Eq. ( |3.5| ) reduces to 



L - ^ {□/[h + d\dJ [uX] - d v f w )} = T [H , (3.13) 

which is still invariant under the transformation (|3.10|) [13]. Thus, we can impose 
the condition 

d v fW = o. (3.14) 

Finally, under the conditions ( |3.11| ), ( |3.12| ) and ( |3.14j) , the field equations of f^ u ) 
and become 

(3.15) 



-kTi 



□/(/«/) 



(3.16) 



where we have defined 1/A == —c\ + (4/9)c3 7^ 0. From Eqs. (|3.11| ) and 03.15 ), we 
find that the symmetric part of T^ v satisfies the conservation law 



d v T, 



(jw) 



0. 



(3.17) 
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and the antisymmetric part of T M „ satisfies 

d v T [H = 0, (3.18) 



which follows from Eqs. (|3.14 ) and ( ft-16 ) |13| 



Let us consider the plane wave solutions of Eqs. ( |3.15|) and (|3.16|) with = 



f { ^{x,x°) = U^ kx +U^ kx , (3.19) 
/ [H (*, x°) = V^e ik ' x + % u e~ ik ' x , (3.20) 

where k ■ >M . Here, and V^ v are constant amplitudes, U and V are their 

complex conjugates, and k^ is a constant wave vector, which satisfy the relations 

W = 0, U^k» = 0, V^r = 0. (3.21) 

Following the prescription given in Section 35.4 of Ref. we impose the transverse- 
traceless gauge condition 

U^C = 0, W„ = 0, V^C = 0, (3.22) 

where is a constant time-like vector. We see that the number of physically 
significant components of f^ u ) is two, while that of is one. 

We next calculate the energy-momentum of the plane waves given by Eqs. ( |3.19 ) 



and (|3.20|) . The dynamical energy-momentum density G T t M of gravitational field 
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has, to lowest order in /„„, the expression 



Grp p 



,(o)a» 



-d°f^dJ (Xp) + dFfWdJw + \d°fdJ+2d°f^d p f [Xa] 



+ d x f^d p f [Xa] --^f^d a f [Xp] + 2d x f^d p f [Xa] 



2e 



(0)1/ 



+ ^dJ^dJ-dPfWdrfw ~ d»f^d a f L 



pv] 



_^ f W\ d J {pv) _ d , f [pa) d J (pu) _ dPf [«A d J {pv) 



- 2d°f^d v f [pa] - d*fMd a fw + d p f^dJ L 

- d p f^d a f [pu] 



(3.23) 



By using Eqs. (|3~T9|)-(|3~2l), G T t M is found to be given by 



Grwi p. 



-2e 



(0)1/ 



2k 



(U pa U pa e 2ikx - 2U pa U pa + U pa U pa e- 2lk x ) 



k»k v 
+ ^2\ 



{V pu V, 



,2ik-. 



pa<- 



2V pa V pa + V^e 



-2ik-x 



pa^ 



(3.24) 



Taking the average of G TJ M over a space-time region much larger than \k\ 1 , we 
obtain 



/G/7-1 M 



,(0)!/ 



r A Up h, UPU 

^(|W n | 2 +|W 12 | 2 )+4^|V 12 | 



(3.25) 



where we have chosen the direction of the space components k of the four vector k p 
as the third axis. The term in the square brackets of the right-hand side (r.h.s.) of 
Eq. ( |3.25| ) is identical to the corresponding term of the canonical energy-momentum 
density t p defined by Eq. (|A.7|) in GR if the antisymmetric part V12 is vanishing.^ 
15 See, for instance, Section 10.3 of Ref. fl5| 
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4 Quadrupole radiation from point masses 



The retarded solutions of Eqs. ( p. 151 ) an d (|3.16|) are given by 



T ( w)(x / ,a: — |sc — aj'|) 



w». *°) = ^ y ^ ^ — - . (4.2) 

respectively. We consider a system of the Newtonian point masses {fn a ,^} (a = 
1,2,..., N) as a gravitational wave source, where m a and denote the mass and 
coordinate of the ath point mass, respectively. For this case, the antisymmetric part 
of the energy-momentum density becomes 

T [H (x,x°) = 0. (4.3) 



Thus, vanishes, i.e., /[^ = 0. Applying the retarded expansion to Eq. ( flTD 
and using the conservation law ( 3.17]) , we obtain the quadrupole radiation formula 



in the rest frame of the system, 

f {af3) (x, x°) = £^dod J d 3 x'x ,a x^T m (x', x°-r), (4.4a) 

f {0a) (x,x°) = -yf {a p)(x',x°-r) } (4.4b) 

N a ^ 

f m {x,X°) = -^^ m aC 2 + -^f(a(3){x',X° - r), (4.4c) 



with r = (x a x Q )2.0 Since we regard each velocity of the mass points to be much 
smaller than the velocity of light, we can use the quantity 



N 



T (00 ) (x, x°) = m * c26i3) ( x ~ *«(*)) ( 4 - 5 ) 



0=1 



16 In our convention, letters from the beginning of the Greek alphabet, a, /3, 7, • • • , and those 
from the beginning of the Latin alphabet, a, b, c, • • • , take the values 1, 2, and 3, unless otherwise 
stated. Also, here we have used the usual summation convention for repeated indices. 
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as the energy density of the system, with t = x jc. Substituting Eq. (|4.5| ) into 
Eq. ( pi~4D , we obtain 

/(a/3) = z — ' /(o«) = Z — ' (4.6a) 

87rr r 87rr 

N a 

III X ™ ry JU Ju t\j 



Anr 87rr 

a=l 

Here, D a| g is the quadrupole moment of the mass distribution defined by 

N 

D a0 = J2 m ^ a , (4-7) 



a=l 



where we have defined D a p == d 2 D a p/ d 2 t. The wave form given by Eqs. ( |4.6j ) and 
( |4.7| ) is the same as the corresponding wave form in GR, which is consistent with 
the result given in Ref. [|Il[ . 



5 Emission rates of energy- momentum and angu- 
lar momentum 

In this section, we examine time averages of emission rates of the energy-momentum 
and angular momentum carried by the quadrupole radiation given by Eq. fl4.6|) . 
In ENGR, for the asymptotic conditions ( [2.45 ) and ( ^.46 ), we have four quantities 



that are conserved as long as they are finite, the dynamical energy-momentum 
"spin" angular momentum Ski, canonical energy- momentum M c „, and "extended 
orbital angular momentum" L " 0. The dynamical energy momentum does not 
depend on the asymptotic form of the Higgs-type field ip k explicitly, but the quantities 
Ski,M c ^ and L " depend on the asymptotic form. With this in mind, following 
Ref. ||, we examine the case in which the asymptotic form of the Higgs-type field 
is 

i) k ~ e^ )k ^ + V (0)fc , 
with constant. We also consider the slightly generalized case in which 

^ k ~pe^\x» + ^ k . 
The asymptotic form of the function p is determined in §BT2 . 
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5.1 The case i/j k ~ e {0)k ^ + ip^ k 
5.1.1 Dynamical energy-momentum loss 

In order to evaluate the emission rate of the dynamical energy-momentum, we inte- 
grate the differential conservation law (|2.37|) over a large solid sphere V with radius 
r, yielding 

d [ tot T k °d 3 x = - [ tot T k a r 2 n a dn , (5.1) 



Jv Js 

where S and dQ represent the two-dimensional surface of V and the differential solid 
angle, respectively. Also, n a stands for the unit radial vector defined by n a = f x a /r. 
Taking into account the fact that the energy-momentum density of point masses 
vanishes for very large r, we can rewrite the r.h.s. of Eq. ( |5.1|) as 

- f tot T k a r 2 n a dQ = - [ G T k a r 2 n a dQ . (5.2) 
Js Js 

The density G T k ^ takes, up to terms of order 0(l/r 2 ), the form 

= e (°>G T (0U t (5.3) 

with 



Or 



(5.4) 



Using the solution ( |4.6|) , and averaging over one period of motion of the system of 
point masses, we obtain 

dE\ G 1 1 \ 

1£J ~ h& \ D <# D <* ~ s DaaD ^/ ' (5 - 5) 

for a total energy E = f —M(q\ of the system, where (■ • • ) denotes the operation of 
averaging over one period of motion of the system of the point masses, and we have 
set e^ k ^ = 5 fe for simplicity. Also, we obtain 

dMa n ^ a\ 

— = 0. (5.6) 

It is worth mentioning that the quantity G T^°\ U is identical to the r.h.s. of Eq. ( |A.14j ), 
and the time average of the emission rate of the dynamical energy-momentum, 
which is given by Eqs. ( |5.5|) and ( |5.6| ), is identical to that of the canonical energy- 
momentum in GR. 
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5.1.2 "Spin" angular momentum loss 

The density G S kl fl has the expression 

where we have defined 



(5.7) 



+ dffMdJw - \d°J^\ v dJ + \&>f {va) drj , (5.9) 

kZ (2W def ^AJM _ ^/(-A)) ^ + _ \grj^ ^ (g 1Q) 

kZ (3)mA CT def 9 AJ( CTM ) _ ^J(aA) _ 

Also, ^ll denotes the energy-momentum density introduced by Landau and Lifshitz, 
whose explicit form is given in Appendix |Aj. Using a method similar to that employed 
in § |5.1.1| , we obtain 

d [ tot S kl °d 3 x = - [ G S kl a r 2 n a dtt . (5.12) 
Jv Js 

In order to estimate the r.h.s. of Eq. (|5.12 ), we set 



^ = + V (0)/1 + ^ , (5.13) 
where ?p( ^ is a constant. From Eqs. Q4.6|) and (|5.7|) - (|5.13|) , we can show thatQ 



dS (P)a \ __ J,(0)a G /- ~ 1 



-^,») (5,4) 



17 Note that we have distinguished between Latin and Greek indices here. (See footnote || on 
page §.) 
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if 

lim = 0, lim ^ = 0, lim rfr g = 0, (5.15) 
where again we have set e'' ^ = 5 k fl . Also, we have 

(/ ' S '"' A ^A^Din-KDav) (5.16) 



dt I 5c 5 

if 



lim 2p a = 0, lim = 0, lim n/^ = 0. (5.17) 



r — >oo r — >oo r — >oo 



The time average of the emission rate of the "spin" angular momentum S a b is the 
same as that of the space-space component of the angular momentum in GR. 

5.1.3 Canonical energy- momentum and orbital angular momentum 
losses 

In the weak-field approximation, the canonical energy- momentum density of the 
gravitational field G T fl u becomes 

G T; = V P ,Z^r M - G T { \ V ^\ . (5.18) 
Then, using Eqs. Q4.6|) , (|5.11|) and ( |5.13| ), we can show that 



dt 

if the conditions 



^=0 (5.19) 



lim fr v = 0, lim < 00 = 0, lim Bv = (5.20) 

r— «x> ' r— >oo ' r— >oo 

are satisfied. 

Finally, we examine the "extended orbital angular momentum." In the weak-field 
approximation, the "extended orbital angular momentum" density of the gravita- 
tional field G M,, uX becomes 



G 



M ux = 2Z (4)^x _ 2Z^ x ^ Vap ^ - 2x u Z^ Xr °i lpa V^ + 2x u Z^ a x r, p , (5.21) 
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where we have defined 



+ 2T 



+ (d p f (aX) - d x f^)) ( d J (ap) - L p dJ 



(5.22) 



4 

+ d x f^dJ (pT) - ~d x fdJ+d T f^d T f (pa) - d T j^dJ [pT) 

- d x f^d T f {pa) - ~d T ff> x \ a d T f+ l -d x f {uT) d T f . (5.23) 

As can be shown by using Eqs. (pj), (^T3|) , (|5T2l| ), ( gggp and ( ggg ), we have the 
relations 



^ = n ^L = n 
<iL 18 \ G 1 ■■ 



(5.24) 
(5.25) 



if the conditions 



lim r^ M „ = 0, lim rip^ 00 = 0, lim r 2 ^ ~ = (5.26) 

r— >oo ' r— >oo ' r— >oo 

are satisfied. The emission rate of the component L£ is finite. However, for the an- 
tisymmetric part of L a p, which is the three-dimensional orbital angular momentum, 
we have 



dL[ a p] 
dt 



0. 



(5.27) 



To summarize, the emission rates dM c fl /dt and dL^/dt are both vanishing if 
the condition ( |5.26| ) is satisfied.^ 



5.2 The case t/j k ~ pe 



(0)k 



+ ^ 



Also in this case, the expressions ( |5.5| ) and ( |5.6| ) for the time average of the emission 
rate of the dynamical energy-momentum hold, as seen from Eq. 



18 Note that £[^1 corresponds to the generator of the Lorentz coordinate transformations and 
that it possibly represents the four-dimensional orbital angular momentum. 
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We express ^ as 

^ = p(r,t)x M + ^ (0)M + ^ . (5.28) 

Then, the expression ( 5.14 ) for the time average of the emission rate of the time- 
space component of the "spin" angular momentum holds if the condition ( 5.15| ) is 



satisfied by ip^ in Eq. (|5.28|) . For the space-space component S a b, using Eqs. fl4~6| 
and (|57?|)- (^T2D , we find 

'dS ab \ 2G3 + p, 



dt 



5c 5 4 
if the condition Q5.17| ) and 

lim p = p c , 

with p c constant, are satisfied by p and ip^ 1 in Eq. G5.28|) . 

From Eqs. (pD, Q, (PH) , (|5~T8|) and ( ggp , we have the relations 

'dM c , 



(5.29) 



(5.30) 



dt 

dM c t 
~dt 



G . 1 1 







if the condition ( |5.20| ) and 

lim p = p c , lim rp = 0, lim rp ; oo = 



(5.31) 
(5.32) 

(5.33) 



are satisfied. 

Under the condition ( |5.30| ), the time average of the emission rate of the "extended 
orbital angular momentum" diverges. However, using Eqs. ( ^4.6| ), ( |5.21| ), ( |5.22| ), 
(|5.23| ) and (|5.28|) , we can show that 



dL 



[Oa] 



dt 

dL[ aj3 ] 

dt 



0. 

2G 1 - p c 
5c 5 4 



\ryLJa-y 



(5.34) 
(5.35) 



if the conditions ( |5.26|) and ( |5.30|) are satisfied. 

For a space-time satisfying the asymptotic conditions ( |2.45| ) and ( |2.46| ), the sum 



Ski + e 



(0)M 



1H 



is well-defined and conserved for the case ip k 



Pce^i+x 11 + 
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^(o)fc ^ ^ 1), as described in Ref. 0. In the case under consideration, we have 



d ,„ , A _ . (o) / dM m 



-< 5 <s m .+w; = -* w ^^ a ;. (5.3 6 ) 

(S'aft + ^[a6])^ — 7^5 (pa-t'Db-j ~ D bl Dajj (5.37) 

if the conditions 

lim = 0, lim r^ v = 0, lim r^ 00 = 0, lim r 2 fr Bv = (5.38) 
are satisfied. Here, again, we have set e^^ k = 5^ k . 

6 Summary and discussion 

In an extended, new form of general relativity, we have examined energy-momentum 
and angular momentum carried by gravitational wave radiated from a system of 
Newtonian point masses in a weak-field approximation. The results are summarized 
as follows. 

1. The form of the gravitational wave is identical to the corresponding wave form 
in GR, which is consistent with the result in Ref. |TT||. 



2. The average value of ( G 7] M ) for the dynamical energy-momentum of a plane 
wave is obtained from Eq. Q3.25 ), and it is the same as that of the corresponding 



canonical energy- momentum in GR fll5|. 



3. The dynamical energy- momentum density G T k ^ takes, up to order 0(l/r 2 ), 
the form given by Eq. ( |5.3| ) with Eq. ( |5.4j ), and this is essentially the same 
as the corresponding expression ( |A.14|) for the canonical energy-momentum 
density in GR, and the time average of the energy-momentum emission rate 
is given by Eqs. ( |5.5|) and (|5.6|), which is identical to that in GR. 

4. The time average of the emission rate of the "spin" angular momentum is 
given by Eqs. (|5.14|) and ( |5.16|) if the conditions ( |5.15| ) and ( |5 . 1 7|) for the form 



( |5.13| ) of the Higgs-type field are satisfied. The expression ( |5.16| ) is the same 



as the corresponding expression for the angular momentum in GR. 



24 



5. The emission rates of both the canonical energy- momentum and the orbital 
angular momentum vanish if the conditions (]5.20 ) and ( |5.26| ) for ip k given by 



the expression (|5.13|) are both satisfied. 

6. Under the condition ( |5.30|) for the form ( |5.28| ) of the Higgs-type field, the 
time average of the emission rates of the "spin" angular momentum and 
the canonical energy-momentum depend on the constant p c . Moreover, the 
time average of the emission rate of the "extended orbital angular momen- 
tum" diverges. However, the time average of the emission rate of the sum || 
Ski + e *; e z-^[H * s finhe, and its space-space component is identical to the 
corresponding expression for the angular momentum in GR if the condition 
Q5.38D is satisfied. 

Finally, we would like to add the following: 

A. As we have stated repeatedly, the dynamical energy-momentum and "spin" an- 
gular momentum are generators of internal translations and internal SL(2, C) 
transformations. The former does not depend on the non-dynamical field ip h 
explicitly, but the latter does. For vierbeins possessing asymptotic forms sat- 
isfying Eqs. (|2.45|) and (|2.46|) , they give the total energy-momentum and total 
(=spin+ orbital) angular momentum of the system when the field ip k is cho- 
sen as ip k — e mx' 1 + ip^ k + 0(l/r^). The generator of the affine coordinate 
transformations, contrastingly, vanishes. The results summarized in above 
are consistent with this. The discussion in Appendix [B| gives further support 
to the choice of the asymptotic form of ip k given by Eq. ( |5.13| ) satisfying the 



condition (|5.17| ). 

B. The "spin" angular momentum depends on the Higgs-type field ip k , and it is 
meaningful if ip k satisfies a suitable condition requiring this field to be the same 
as the Minkowskian coordinates on the boundary of a sphere of infinite radius, 
as is known from §[| and from Refs. || |7|, || [J. The field if) behaves like a 
Minkowskian coordinate system under internal Poincare transformations, and 
its existence is a necessary consequence of the structure of the group P and a 
basic postulate regarding the space-time. Also, this field is closely related to 
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the existence of the spinor structure Q . However, the physical and geometrical 
meaning of this field has not yet been fully clarified. 



C. In considering the energy-momentum and angular momentum, there are two 
possibilities in choosing the set of independent field variables 0, |6|, i.e., the set 
{ip k , A k ^, <f> A } and the set {ip k , e fc M , 4> A }. In the present paper, we have employed 
{ip k , A k ^, cj) A } as the set of independent field variables, because this choice is 
superior to the other, as shown Refs. || and || .0 



D. (a) The asymptotic conditions (|5.26|) and (|5.30| ) for the field ip k are stronger 



than the corresponding conditions for vierbeins whose asymptotic forms 
satisfy Eqs. ( j2.45| ) and ( |2.46| ) 0. This is natural, because vierbeins in 
which gravitational waves propagate do not satisfy the conditions ( 2.45 ) 
and ( |2.46| ) in general. 



(b) For the form ( |5.28| ) of a Higgs-type field, the time average of the emission 
rate of the canonical energy, which is given by Eq. ( [5.31| ), is identical to 
the corresponding expression in GR if the condition ( |5.33j ) with p c = is 
satisfied. However, when p c = 0, neither the time average of the emission 
rate for the space-space component of the "spin" angular momentum 
nor that of the orbital angular momentum are the same as the corre- 
sponding expression in GR. The time average of the emission rate for the 
space-space component of the sum Sf-i + ^e^ is the same as the 

corresponding expression in GR. However, the sum Sf-i + e^ >tl k e^ lJ l Ly^ v -\ 
is an artificial quantity, because the "spin" angular momentum Ski and 
orbital angular momentum Li^j are associated with transformations that 
differ with each other. 

E. The transformation property of the dynamical energy-momentum density of 
gravity in ENGR is different from that of the canonical energy-momentum 
density of gravity in GR. The former behaves as a vector density under general 



19 It should be noted that when {ip k , A k ^, (f> A } is employed as the set of independent field variables, 
the Lagrangian L is considered to have an explicit if) dependence, because L is a function of 
tp k „ + A k <j) A and their derivatives. 



p k _ „/,fc 
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coordinate transformations, while the latter is not tensorial. Nevertheless, the 
two densities take the same form up to order 0(l/r 2 ). 



F. In ENGR, the world line of a macroscopic test body is the geodesic of the 
metric g, as in the case of GRj^| and hence ENGR with the condition ( |3.12| ) 
accurately describes the variation of the period of motion of the binary pulsar 
PSR1913+16 [IE]. This is consistent with results in Refs. [IU| and [O]. 



G. As far as the gravitational radiation from Newtonian point masses treated 
in the weak-field approximation is concerned, the losses of energy-momentum 
and angular momentum as well as the wave form are independent of the pa- 
rameter C3 in the gravitational Lagrangian density, and they are identical to 
the corresponding quantities in GR. The effects of C3 reveal themselves at 
higher orders. 



H. We have considered the case of a weak field under the condition ( 3.12 ), but 
the field equations ( |3.4| ) and ( |3.5| ) can be solved also in the case with the 
following Jr7| : 

Cl + c 2 ^0, c x --C3^0. (6.1) 
y 



The dynamical structure of the system with the condition ( |6.1| ) is significantly 
different from that with the condition (|3.12|) . Although the values of the 
parameters c\ and c 2 are severely restricted as 

ci ~ -— ~ -c 2 (6.2) 

6K 



by the results of Solar System experiments [O, there still remains the possi- 



bility that Eq. (|3.12|) is not satisfied. The gravitational radiation for the case 
in which the parameters Cj, c 2 and C3 satisfy Eq. (|6.1|) is also worth examining. 



From [A], |C| and [D], we deduce the following: The choice {ip k , A k , <ft A } with 
Eq. ( |5.13| ) satisfying the conditions ( p . 1 5| ) and ( p.!7| ) is superior to the other possible 



choices, and the generator of internal Poincare transformations accurately describe 



2(1 



Note that the behavior of the Dirac particle in ENGR differs from that in GR. (See Ref. fL3|.) 



27 



the energy-momentum and angular momentum for a wide class of gravitating sys- 
tems, including space-times in which there are gravitational waves. 



20, 21 



In the teleparallel theory of gravity, there have been several attempts [Tj| [19, 
73|1 to define well-behaved energy-momentum and angular momentum 



densities. For the case of the teleparallel equivalent of general relativity, i.e. the case 
with the condition Q2.27p in our notation, this problem is studied in Refs. 



191 BO 



The gravitational energy-momentum density hj a p in Ref. is the same as our 



G 



± k 



. In Refs. |19| and p0| , Hamiltonian formalism is developed, and a natural 
definition of the energy-momentum density of the gravitational field is given. In 
addition, the angular-momentum density is examined in Ref. IT 



In Ref. 21 



an 



energy-momentum current that transforms as a tensor under diffeomorphisms of the 
space-time manifold and under global SO(l, 3) transformations is proposed in a co- 



frame field formulation of the general teleparallel theory of gravity. In Refs. |22] and 
|23| , the energy-momentum and angular momentum densities for general teleparallel 
theory of an isolated gravitating system are examined. 

The energy-momentum and angular momentum for gravitational waves, however, 
are not examined in Refs. f[ 



lj, |fj, 0, |2[ H, and in Refs. [[TJ, |fj, |T|, ||, [23 



the energy-momentum and angular momentum are not related to the generators 
of internal Poincare transformations, and there appears no field that corresponds 
to our ip h . It is worth examining the relation between their [19, [Hj, 
energy-momentum and angular momentum densities and ours. 
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A Linearized Einstein Theory 



For convenience, we give here a short summary of linearized GR. 
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The gravitational Lagrangian density is given by 



r def 1 r 
L G R = 



y. p) [u A J [p v) [X p 



where we have defined the Christoffel symbols by 



pu\ ^ \ gXP ( 9 » 9vp + 9v9pfl ~ dp9 ^ ' 



The Einstein equation takes the form 



(A.l) 



(A.2) 



Rfiu({}) - ^9nvR{{}) = kT^ , 



(A.3) 



where we have defined the Riemann- Christoffel curvature, Ricci tensor and scalar 
curvature by 



X p) [a is ) [X v) [a p 



(A.4) 

(A.5) 
(A.6) 



respectively. Also, denotes the energy-momentum density of the gravitational 
source. The canonical energy- momentum density is defined by 

dL GK 



def 



t„ = 6„ L GR - 



9pu, P 



(A.7) 



We consider a metric perturbation h^ v from Minkowskian space-time, i.e., 



(A.8) 



The linearized field equations are given by the form in Eq. ( [3.15D , up to an overall 
factor of 2, 

-2kT^ , (A.9) 



Dh 



with the harmonic coordinate condition 



d v hT = 0, 



(A.10) 



29 



where we have defined 



i>;,u = - -rj^uh , h = rf v h 



(A.ll) 
(A.12) 



Note that the perturbation h^ u corresponds to Zf^u) in §H At lowest order, t " 
becomes 



■iiit;; = <\;' ( \d p h XtT d x h po - -o x h>"<) x h p ,, + UrinU, 



which is approximated, up to order 0(l/r 2 ), as 



(A.13) 



(A.14) 



Finally, the energy-momentum density of gravity proposed by Landau and Lif- 
shitz |14| is defined by 

(-gKl = 9^ - (-g)T^ , (A.15) 

where 9 pu is the symmetric energy- momentum density defined by Eq. (j2.51|) . By 
using Eqs. (|2.51|) and ( |A.3| ), we obtain from Eq. ( |A.15| ) the following: 



Ik t hh 



a 



T 



+ g vX g pa 



+ g xp g aT ({ 



A pj [a T) [X T J [per) [X o) ip T 

V 1 f T 1 f V 1 f T 1 f^llT 



A r J \p a j [p a j [X rj ^ct rj [A pj [A pj [a r 



v 



A* - + - 
A rj ctJ [p tr J [A r 

A aj [p t) \X p) [a r 



p 

A pj l^a 

(A.16) 



To lowest order, this takes the form 

2K{-g)t& = rf v (\d x h^d p K x + l&>hdji - \d"h aX d p h aX \ + x -d^d v h pa 



l -^hd v h - d»h pa d p h UCT - d p h^d v h pa + d p h pa d p h v 



(A.17) 
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B Angular Momentum Loss Derived from 
Energy- Momentum Loss 



Following the argument given in problem 3 in Section 110 of Ref. ||14|| , we derive 
the time average of the orbital angular momentum loss for Newtonian point masses 
from the time average of the dynamical energy loss given by Eq. (|5.5|) . The result 
supports the discussion of the "spin" angular momentum loss given in § |5.1.2| . 

We represent the time average of the energy loss of the system as the work of 
the "frictional forces" q acting on the Newtonian point masses: 



dE* Y 



dt 



£(<?«.•*.). (b.i; 



a=l 



The time average of the loss of angular momentum, I = f ^2^=1 £a X f^aia, is given 
by 

( " x - £ ((€. x Qa ) a ) = £ eof* (&Ql) , (B.2) 

a=l a=l 



dt 



where the symbol e Q( g 7 denotes the three-dimensional antisymmetric tensor with 
ei23 = e 123 = 1- Note that a Newtonian point mass in ENGR is subject to Newton's 
equation of motion. To determine g a , we write Eq. fl5.5|) as 

■fiT/ = ~-^E ('D a pD a f3 - -DaaDpp) = (b^D^ - -D aa D { ^\ , (B.3) 



where we have used the fact that the average values of the total time derivatives 

a/3 



vanish. Here, D^l represents the fifth-order derivative with respect to t. Substitut- 



ing D a p = J2 a =i m aH2ia + into Ec L- Q) and comparing with Eq. ( pT|) , we 

find 

Substitution of Eq. flB.4|) into Eq. ( |B.2|) gives the result 



= -ge a/37 ^ _ lsy D lJ) = -g^ (i),^) • (B.5) 

This is equivalent to Eq. ( |5.16| ), with Eq. ( |5.13| ) satisfying the condition ( p.!7| ). 
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